Doping dependence of magnetic excitations of ID cuprates as probed by Resonant 

Inelastic x-ray Scattering 
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We study the dynamical, momentum dependent two- and four-spin response functions in doped 
and undoped ID cuprates, as probed by resonant inelastic x-ray scattering, using an exact numerical 
diagonalization procedure. In the undoped t — J system the four-spin response vanishes at tt, 
whereas the two-spin correlator is peaked around 7r/2, with generally larger spectral weight. Upon 
doping spectra tend to soften and broaden, with a transfer of spectral weight towards higher energy. 
However, the total spectral weight and average peak position of either response are only weakly 
affected by doping up to a concentration of 1/8. Only the two-spin response at n changes strongly, 
with a large reduction of spectral weight and enhancement of excitation energy. At other momenta 
the higher-energy, generic features of the magnetic response are robust against doping. It signals 
the presence of strong short-range antiferromagnetic correlations, even after doping mobile holes 
into the system. We expect this to hold also in higher dimensions. 

PACS numbers: 78.70.-g 74.72.-h 78.70.Ck 71.27.+a 



I. INTRODUCTION 



Cuprate materials have proven a fertile ground for the 
study of strong electronic correlations and quantum mag- 
netism. Spectroscopic techniques are powerful tools to 
obtain information about these properties, since they 
provide direct information on the electronic and mag- 
netic elementary excitations, which are related to the en- 
ergy spectra, the crystal structure, and so on. Among 
the spectroscopic techniques used to probe magnetic ex- 
citations in cuprates, Resonant Inelastic X-ray Scatter- 
ing (RIXS) has gained much interest because of the re- 
cent increases in energy and momentum resolution, due 
to the enhanced brilliance of synchrotron x-ray sources 
and the advances in instrumentations. This has placed 
RIXS at the forefront in the study of the momentum de- 
pendent electronic and magnetic responses over a wide 



energy range 
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In the RIXS process, x-ray radiation is inelastically 
scattered by the matter and the change in energy, mo- 
mentum and polarization can be related to intrinsic exci- 
tations in the material. This process is resonant because 
the energy of an incoming photon is tuned to match an 
element absorption edge, thus allowing a large enhance- 
ment of the scattered intensity. In this way the x-ray 
photons can couple to charge, spin and orbital degrees of 
freedomJ^^ 

In the context of the cuprates, it is by now well- 
established that RIXS can detect the momentum depen- 
dence of charge excitations that are related to the elec- 
trons and holes in the d shell,— i 24 i 25 , but it has also been 
proved, both experimentally and theoretically, that RIXS 
is sensitive to the magnetic excitations of cuprates. 



A. Magnetic RIXS in cuprates 

In the magnetic sector, RIXS can both create single 
and double spin flip excitations, corresponding to singlc- 
and bi-magnon excitations in ordered Hcisenbcrg antifcr- 
romagnets (AFM). Moreover, it has been also predicted 
that threc-magnon scatterings contribute substantially to 
the magnetic spectral weight,— whereas singlc-magnon 
excitations were observed in direct RIXS experiments at 
the Cu L 3 edge on a thin film of I^CuC^r ^ 21 ' 28 and 
also on small crystals of S^CuC^C^.— 

From a theoretical point of view, when the magnetic 
moment lies in the plane of the x 1 — y 2 orbital, direct 
spin-flip scattering in cuprates is allowed for symmetry 
reasons. Thus, at least for this class of materials, L-edge 
RIXS can be placed on the same footing as neutron scat- 
tering, and both are related to the two-spin dynamical 
correlation function^ 

At the Cu K-edge, the RIXS process is indirect and 
single-magnon scattering is forbidden. In this case, 
magnetic excitations turn out to be due to bi-magnon, 
as observed in insulating and doped L^-zSr^CuC^ 
and Nd2Cu04i 16,22 Subsequently, excitations with a bi- 
magnon-likc dispersion have been observed in cuprates 
with high- resolution L-edge RIXS,— ^ and M-edge 
RIXSi^ This makes RIXS thus complementary to optical 
Raman scattering, which also measures the bi-magnon, 
but only at zero momentum transfer £2r£L. 

The microscopic mechanism by which the bi-magnons 
couple to the intermediate state core-hole in Cu K- 
edge RIXS is by the core-hole locally modifying the su- 
perexchange constant. This leads to the measurement 
of a four-spin correlation function that can be derived 
in detail via the Ultrashort Core-hole Lifetime (UCL) 
cxpansion i 26 i 27 ' 38 i 39 This approach yields a momentum 
dependence of the cross section in agreement with ex- 
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pcrimcnts on undopcd cupratcs in the Neel state. In 
particular it reproduces the lack of intensity at g=(0,0) 
and (-7T, 7r) . 16 i 22 At the Cu L- and M-edges, the coupling 
mechanism to bi-magnon excitations is similar, resulting 
in the same cross section to lowest order in the perturb- 
ing effective potential between the spin and the core-hole, 
due to the fact that the superexchange in the intermedi- 
ate state is, in this case, not just altered but completely 
blocked locally^ 

Despite the success in describing undopcd cuprates, 
the doping dependence of both single- and multi-magnon 
measured by RIXS is rather little unexplored, however 
with the first data available showing intriguing behav- 
ior. For K-edge measurements on L^-xSr^CuCj, the 
intensity of the spectral features is found to generally de- 
crease as x increases . 16 ' 22 The persistence of the highest 
peak at (n,Q), located around 500 meV, for doping val- 
ues x = 0.07, well into the superconducting phase, shows 
that this excitation survives even if long-range magnetic 
order is absent, as long as significant short-range mag- 
netic correlations are present, which is well-known to be 
the case in the superconducting state of LSCOi 40 ' 41 The 
evolution upon doping observed in L-cdgc RIXS mea- 
surements on LSCO is even more captivating, showing 
the existence of a high energy "undopcd" branch in ad- 
dition to the lower-energy dispersive features measured in 
neutron scattering, which can be related to the presence 
of a stripe liquid 21 -. 



B. Aim and Outline 

In this paper we analyze the doping dependence of 
magnetic RIXS in cupratcs, and particularly we explore 
the evolution of the magnetic excitations with doping. 
We focus on ID systems because i) strong quantum fluc- 
tuations are present, that maximally effect the magnetic 
ordering, ii) the existence of an exact theoretical results 
for magnetic excitations in the undopcd ID Hciscnberg 
AFM, formulated in the spinon language, provides a 
stringent reference test for our numerical spectra, iii) the 
sampling of the Brillouin Zone (BZ) is much more dense 
than in the 2D case, and iv) the results are directly rel- 
evant for RIXS on ID cuprates, such as Sr 2 Cu03, where 
experiments probing the low-energy magnetic excitations 
are entirely feasible^ These ID spin liquid states may be 
relevant to the proposed stripe liquid behavior in high- 
temperature superconductors. 

We present numerical calculations of different-time 
Two-Spin (TS) and Four-Spin (FS) correlation functions, 
evaluated on 22 sites chain for the Hciscnberg model and 
16 sites chain for the t — J model, in the low doping 
regime with a doping concentration up to 1/8. Even if 
the analysis is limited by the finite length of the system, 
one can extract relevant information about the intensity 
and the dispersion of the magnetic features by quantita- 
tive analysis of total spectral weight (Wo), and its first 
moment (Wi), corresponding to the average energy peak 



position. Besides, the higher energy features are expected 
to be less affected by finite size effects. 

We firstly consider the undoped case assuming a de- 
scription based on the nearest-neighbour AFM Hciscn- 
berg model on a chain. It turns out that both TS and FS 
correlation functions detect two-spinon excitations, with 
a spectral weight that is concentrated at the lower bound- 
ary of the continuum of excitations. We point out that 
TS and FS access excitations belonging to orthogonal 
subspaces, having total spin S=l and S=0, respectively. 
Moreover, we recover crucial differences about momen- 
tum dispersion and selection rules: TS is peaked at tt 
where FS vanishes, whereas the latter is peaked at ir/2. 
This result highlights the different length-scale between 
the TS and the FS excitation. The latter has a charac- 
teristic length-scale of 2a, where a is the lattice spacing, 
as two (neighboring) exchange bonds are broken in the 
intermediate state. This shifts the maximum from tt to- 
wards ±7r/2. These results arc consistent with recent 
Bethe Ansatz calculations^. We subsequently study the 
spectral evolution upon doping, by using the same ap- 
proach on a t — J chain containing mobile holes. The 
magnetic response is found to soften and broaden as a 
function of the doping concentration. Nevertheless, both 
the momentum dispersion and the total spectral weight 
are only slightly affected by a doping of 1/16 and 1/8. 
Only the two-spin response at tt changes strongly, with 
a large reduction of its spectral weight and an enhance- 
ment of the excitation energy. At other momenta the 
higher-energy, generic features of the magnetic response 
are robust against doping. This shows that the strong 
short-range antifcrromagnetic correlations that are still 
present after doping give rise to higher energy, damped 
magnetic excitations with considerable spectral weight. 



II. MAGNETIC EXCITATIONS IN THE ID t — J 
MODEL 

The t — J Hamiltonian is one of the most studied model 
Hamiltonians in the context of high temperature super- 
conductivity in doped, quasi-2D cupratcs. Naturally, it 
is directly relevant to quasi-lD cuprates as well, among 
which edge-sharing ID spin chains in Li2Cu02 and 
GeCuC>3 and corner sharing in SrCu02 and Sr2Cu03. 
Charge excitations have been extensively studied by ex- 
act diagonalization of Hubbard or charge-transfer mod- 
els relevant for these systems, also in the context of 
RIXS 44 ' 45 . In a Hubbard-type model for an undopcd 
chain the presence of singlet excitations of spinon pairs 
was identified, which is due to the presence of doubly 
occupied sites^ 6 -. Very recently the RIXS response of 
the undoped spin-only chain was computed by Bethe 
Ansatz, which will serve as a benchmark for the present 
exact diagonalization study.—. However, the magnetic 
response of doped t — J chains has received remarkably 
little attention in the literature so far. Studies of spin dy- 
namics in the ID t — J model have considered the limit 
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of quarter- filling, with a hole concentration of 3/4, by 
Monte-Carlo 47 ' 48 , exact-diagonalizatior* 4 ^ and recursion 
methods^, focusing on the regime of weak and strong 
coupling, comparing one- and two-dimensional features, 
as well as determining the evolution of the spinon spec- 
trum in presence of magnetic anisotropy. Here we con- 
sider the weakly doped chain, with doping concentrations 
where in the 2D cuprates superconductivity appears and 
we concentrate on the specific question how far the dif- 
ferent types of magnetic excitations to which RIXS is 
sensitive are affected by the presence of mobile charge 
carriers. 



A. Zero doping Heisenberg limit 

Without doping the ID t — J model reduces to the 
Heisenberg AFM, which is one of the few many-body 
problems where the ground state, which is a SU(2) sin- 
glet, and the lowest excited states are known exactly^ 
It is described by the Heisenberg Hamiltonian 

Hiieia = J Si • Sj , (1) 

(iJ) 

where J is the exchange coupling and Si is the spin op- 
erator at the t-site. It is well-known that this model fails 
to develop long-range Neel ordering (where neighbour- 
ing spins point anti-parallel to each other) even at the 
lowest temperatures but rather it has an algebraically 
decaying spin-spin correlation ! 52 ! 53 Moreover, the basic 
excitations are spinons^^— which are topological exci- 
tations that can be visualized as twists of it in the spin 
order. Spinons are fractional particles that possess spin 
values of S = 1/2, whose dispersion relation is given by 
e(p) = 7r/2|sinp|, p e [— 7r, 0] in unit of J. Because of 
quantum mechanics constraints, only an even number of 
spinons can be created. 

The lowest energy excitations are made of two spinon 
states, that live within a continuum in (k, uj) defined by 
the kinematic constraints of momentum and energy con- 
servation: k = —pi — P2 and oj = e(p±) + e(j)2)- Hence, 
for a fixed external momentum, there exists an interval 
in frequency given by the conditions: 

lu > Ll>2 h{k) = 7r/2| sinfc|;w < U2 u(k) = vr | sinfc/2| , 

(2) 

with k G [0, 27!"]. Here, u>2,L(k) and LU2,u{k) correspond to 
the lower and upper bands, respectively. This behavior 
obviously very different from 2D and 3D systems, where 
the excitations in the ordered state are spin-waves, which 
possess a spin value equal to one and follow a well-defined 
trajectory in energy and wave-vector space in complete 
contrast to the multi-spinon continuum. 

Two-spinon excitations are routinely measured by in- 
elastic neutron scattering 5 - 7 - - — , but the same response 
function can also be measured with RIXS. A wealth of 
theoretical work has been done on the TS correlation 
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FIG. 1. Two-spinon (left panel) and four-spinon (right panel) 
continua, stemmed in the colored regions. Color scale is not 
related to the spectral weight. 



function related to neutron structure factor and the ex- 
act theoretical spectrum perfectly matches with experi- 
mental results £ir— 

The two-spinon part of the TS function is finite inside 
the two-spinon continuum, and by construction it van- 
ishes identically outside of it. It has been demonstrated 
that approximately two thirds of the excitations in this 
band are indeed of the two-spinon type. The remain- 
ing part has a rather small spectral weight and is car- 
ried mostly by four-spinon excitations^ In Fig|TJ these 
two-spinon and four-spinon bands are schematically rep- 
resented. Recently it was proven for the FS response 
function, even if it involves the excitation of two spins, 
overwhelmingly fractionalizes into two-spinon states.— 

Within this context, we now start a systematic com- 
parison of the numerical spectra related to the TS and 
FS correlation functions. We will demonstrate that our 
results recover the main features (e.g., dispersion and se- 
lection rules) of the exact theoretical results, providing 
the starting point for the subsequent discussion about 
doping evolution. 



B. Cross section for magnetic RIXS 

In the Introduction we mentioned that both direct 
and indirect RIXS can probe magnetic excitations. In 
direct RIXS, single spin-flip excitations can be made at 
the 2p — > 3d edges of Cu because of the large spin-orbit 
coupling of the 2p core- hole.— Since in the intermediate 
state spin and orbital angular momentum separately 
are no longer good quantum numbers, orbital and spin 
can be exchanged and direct spin-flip processes can 
in principle be allowed. As in the neutron scattering, 
the RIXS cross section consists of a local structure 
factor (depending on the polarization, the experimental 
geometry and the excitation mechanism) multiplied by 
the appropriate spin susceptibility^ 

S(q,w) oc£|<0|S(<z)|/) 2 <J(a;-w/i) , (3) 
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FIG. 2. A schematic picture for a single spin-flip [S(q, to)] 
and a double spin-flip [0(q, w)] in an AFM Heisenberg chain. 
The first fractionalizes into two-spinon having parallel spins. 
The latter into two-spinons carrying opposite spin (total spin 
S=0) . Spinons are emphasized by the dashed box. The circle 
indicates the sites where the spin flip process occurs. 



where |0) is the ground state, |/) an excited state, w/, 
the energy lost by the photon, and S(q) = J2 i exp igRi S Z i 
is the single-spin form factor. It was shown both 
theoretically ) 19 i 28 and experimentally 2 ^ that spin-flip ex- 
citations are a result of the effect of L Z S Z operator in the 
core hole spin-orbit coupling, so one can get a pure spin- 
flip transition for a Cu 2+ ion with a hole in the 3d x 2_ y 2 
orbital only if the spin is not parallel to the z-axis. 

For indirect RIXS at Cu K-edges (Is — > 4p), the core 
hole couples to the spin degree of freedom locally mod- 
ifying the supcrexchange interactions! 26 ! 27 In this pro- 
cess, the total spin of the valence electrons is conserved, 
and only excitations with at least two spins flipped 
(with total S z = 0) are allowed] 16 ' 22 Detailed calcula- 
tions of the magnetic response functions within the UCL 
expansio n 38 ' 39 demonstrated that the magnetic correla- 
tion function, measured by indirect RIXS, is a four-spin 
correlation one, 

0(q,w)<Kj2mO(q)\f) 2 S(u>-LJ fi ) , (4) 
/ 

where 0(q) = J2i ex P iqRi (J2s^i ' ^i+s) is the two-spin 
form factor. At the transition metal L-cdges, a similar 
mechanism occurs (in addition to the direct single spin- 
flip scattering discussed above) because the photo-excited 
electron in the 3d subshell frustrates the local superex- 
change bonds H 

To calculate S(q, u) and 0(q, w), we apply the Lanczos 
algorithm to a 22-sites Heisenberg chain and we employ 
eigenvectors and eigenvalues for the evaluation of the cor- 
relation functions as in Eqs. In addition, we perform 
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FIG. 3. Inelastic part of S(q, u) at different hole concentration 
associated to the unequal time correlator defined in Eq. [3] 
evaluated for a 22-sites Heisenberg chain and 16-sites t — J 
chain. For the t — J model, the hopping is fixed at t — 3 J. The 
broadening of the lorentzian is 7 = 0.5 J(0.05 J) for the right 
(left) panel, respectively. 
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FIG. 4. Inelastic part of 0(q, uj) at different hole concentra- 
tion associated to the unequal time correlator defined in Eq. 
IH evaluated for a 22-sites Heisenberg chain and 16-sites t — J 
chains. For the t — J model, the hopping is fixed at t = 3J. 
The broadening of the lorentzian is 7 = 0.5J(0.05J) for the 
right (left) panel, respectively. 
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an extensive analysis of the momentum dependence of the 
frequency moments. We evaluate the average peak posi- 
tion with respect to the double differential cross section 



Wi oc / i 



diiduj 



dui and we use the total weight 



dQdu 

Wo oc / %ndH) d*** ^° S am information on the relative ra- 1 ^ 
tio as a function of doping concentration and to compare 
the total intensity of the bare TS and FS spectra. Note 
that these outcomes are possibly easier to compare with 
experiments, with respect to the differential cross-section 
itself. The reason is twofold: i) the presence of statistical _ 
errors in experiments is less relevant for the integrated 
quantities, so a direct comparison with theory is feasible; 
ii) since lineshapes of the theoretical spectra are typically 
not Lorcntzian, the average excitation energy of the cal- 
culated spectra is expected to be the most representative 
theoretical result. 

Before describing the results obtained within our sim- 
ulation, we would like to point out that the magnetic 
excitations described by Eqs. [3JH] belong to orthogonal 
subspaces, as illustrated in FigfJl The figure aims to 
clarify that the S(q, ui) TS response is related to a single 
spin-flip that gives rise to a two spinons carrying parallel 
spin. Instead 0(q,ui) FS response function is due to a 
double spin-flip that fractionalize into two spinons, hav- 
ing opposite spin. Starting from the SU(2) singlet ground 
state, the excitation governing S(q,uj) thus carries S=l 
while for 0(q,uj) the excitation has S=0. This crucial 
aspect allows to distinguish between 'polarized' and 'un- 
polarized' excitations and suggests that the two magnetic 
responses can have a different sensitivity to for example 
spinon-spinon interactions or an external magnetic field. 

In the top panels of Figs. [30 we report the undoped 
inelastic intensity evaluated numerically starting from 
TS and FS correlation functions, respectively. There 
are several common features emerging: the dispersions 
show similar behavior, the excitation spectrum is mainly 
located in the two-spinon band, both for S(q, ui) and 
0(q,u>), with an energy scale u>2,l < uj < u)2,u- More- 
over, the intensity is dominated by the spectral weight at 
the lower threshold of the two-spinon continuum. These 
results confirm that both S(q, ui) and 0{q, ui) mainly frac- 
tionalize in two-spinon excitations*^. From a closer in- 
spection at the energy region comprised between the up- 
per boundary of the two-spinon continuum and the up- 
per boundary of the four-spinon continuum, high energy 
tails emerge as consequence of the four-spinon part of 
the structure factors, that are finite but very small and 
rapidly approaching to zero. 

Nevertheless, spectra do show differences in dispersion. 
Looking at the white circles of Fig. [51 we infer that for 
the correlation function 0(q,u>) the spectrum disperses 
from zero at q = 0, up to a maximum of ~ 2J at 7r/2, the 
intensity is suppressed at ir and nearby, a fact that is due 
to structure factor e lqRi and not to the density of states, 
as reported in lower panel of Fig. [5jwhere the moment W\ 
is plotted. The vanishing of the RIXS intensity described 
by 0(q, ui) at the antifcrromgnctic wavevector 7r is due to 
the cancellation of the sums over the two sublattices in 
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FIG. 5. Total weight (upper panel) and first moment (lower 
panel) of the TS (left side) and FS (right side) correlation 
function on a 16-sites chain, for different doping concentra- 
tions. For convenience of graphical representation, Wo is 
scaled to the total weight at n. Missing points for W\ are 
due to vanishing spectral weights. 



an antifcrromagnctic order. A more detailed explanation 
can be found in Refill. On the contrary, S(q, ui) has its 
maximum at ir (see Fig. [SJ. The appearance of a satellite 
peak at the AFM point located at higher energies can also 
be noticed. Finally, we remark that, as long as the bare 
correlation functions are considered, the total weight is 
larger in the FS case, almost all over the BZ (away from 
7r), while the highest peak intensities are comparable. 



C. Doping dependence 

The doped Hcisenbcrg model is here considered in 
terms of the single-band t — J model^ 



Ht-.r = - t ^2 ( d\, a d l+ s.a + h. 

(i.8),a \ 

+ J^S i -S j , (5) 



(i,3) 



where the sums run over all the bonds, counted 

once, and d operators describe electrons in the d elec- 
tronic levels with the constraint of no double occupancy. 

By following the same procedure described in the pre- 
vious section, we analyze the evolution upon doping of 
the TS and FS spectra. Two cases are considered: the 
hole concentration 71^=1/16 and n^=l/8, respectively. 
In Figs. [3] and 0] we report the spectra for t/J = 3, a ra- 
tio that is typical for corner-sharing cuprate perovskitcs. 
As one can see, the spectra are generally softened and 
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broadened and the intensity of the highest peaks is re- 
duced as n^ is increased. This effect is more dramatic at 
tt for TS, where the peak is strongly damped. In con- 
trast, the doping does not affect substantially the disper- 
sion of the excitation continuum, the only effect being a 
transfer of spectral weight towards the upper threshold of 
the two-spinon continuum, with tails in the four-spinon 
region. Concerning the selection rules, the inelastic in- 
tensity is vanishing at q = for both TS and FS, and it 
is suppressed at tt for 0(q,u). 

Considering the frequency moments reported in Fig. [SJ 
interesting features come out. Namely, as long as the TS 
is concerned, major differences occur around the AFM 
wave- vector, where Wo strongly decreases with the in- 
creases of the doping. This result is somewhat expected 
since AFM correlations arc weakened by doping. The 
evolution of Wi suggests that the average peak position 
stays almost unchanged upon doping; nevertheless it is 
shifted to to ~ 2 J approaching tt. As compared to the TS 
correlation function, the intensity of the FS appears to be 
renormalized in a less dramatic way. The total weight is 
even increased in the first half of the BZ (see Fig. [5]) and 
the intensity of the highest peak at ir/2 is much less soft- 
ened than the corresponding TS at tt. A further look at 
Wi allows to conclude that n/ l =l/16 leaves the dispersion 
of the average peak position unmodified while, by increas- 
ing rih, the dispersion is shifted to slightly higher energies 
(uj ~ 2.5 J). From these outcomes, we can conclude that 
both the magnetic dispersions mapped by S(q, ui) and 
0(q,u>) are quite robust against the doping, as long as 
the lightly doped regime is considered, and we also infer 
that the robustness of the two-spinon feature is directly 
linked to the existence of short-range AFM correlations, 
irrespective of doping concentration. 

Finally, concerning the transfer of spectral weight to- 
wards the upper threshold of the two-spinon continuum, 
we may deduce that it can be related to an 'itinerancy 
effect' due to the fact that the introduction of mobile 
carriers induces charge fluctuations that couple to spin 
excitations. We note that an analogue effect is known to 
be played by correlation energy. For a single-band ID 
Hubbard model at half filling, it has been showed that 
for large values of the on-site repulsion U , the spin cor- 
relations are dominated by virtual hopping processes of 
electrons and are described in terms of a spin- 1/2 Heisen- 
berg chain; in this limit, the spectral weight of spin struc- 
ture factor S(q,uj) is concentrated at the lower spinon 



boundary^ If U is decreased, the real hopping of elec- 
trons becomes important and eventually dominates the 
spin response. The electron itinerancy may influence the 
magnetic correlations, and it is found that the spin struc- 
ture factor of the Hubbard model differs from the spec- 
trum of the Heisenbcrg model since a significant spectral 
weight is concentrated on the upper spinon boundary^ 
We therefore infer that in our case those 'itinerancy cor- 
rections' arising from charge excitations may be driven 
by doping with holes, both on TS and FS functions. 
III. CONCLUSIONS 

For finite Heisenberg and t — J chains we have cal- 
culated the two-spin (TS) and four-spin (FS) correlation 
functions, which determine the magnetic RIXS spectra of 
cupratcs. The calculation is performed by means of the 
Lanczos algorithm applied to 22-sites Heisenberg and 16- 
sites t — J models. We firstly show that both functions 
measure multi-spinon excitations, and we find that the 
dominant contribution arises from the two-spinon con- 
tinuum. Depending on which function is considered, one 
can access to 'polarized' (S = 1) or 'unpolarizcd' (S = 0) 
excitations, resulting also in a different momentum dis- 
persion. Then, as long as the Heisenberg model is con- 
cerned, one observes that O{q,to) peaks at tt/2 while 
S(q, uj) peaks at tt, where 0{q, uj) vanishes; away from tt, 
the first has a total weight lower than the latter. When 
doping is introduced, both the TS and FS spectra are 
softened and spread out. The maximum peak intensities 
are strongly renormalized, but the total integrated in- 
tensity is slightly modified: the spectra tend to broaden. 
The TS and FS spectra show a general transfer of spec- 
tral weight towards the higher energy sectors of the two- 
spinon band. The main features of the TS and FS spec- 
tra, however, exemplified by for instance their zeroth and 
first moment, show a remarkable robustness against to 
doping. 
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